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Abstract 

Context. KIC 4247791 is an eclipsing binary observed by the Kepler satellite mission. 

Aims. We wish to determine the nature of its components and in particular the origin of a shallow dip in its Kepler light curve that 
previous investigations have been unable to explain in a unique way. 

Methods. We analyse newly obtained high-resolution spectra of the star using synthetic spectra based on atmosphere models, derive 
the radial velocities of the stellar components from cross-correlation with a synthetic template, and calculate the orbital solution. We 
use the JKTEBOP program to model the Kepler light curve of KIC 4247791. 

Results. We find KIC 4247791 to be a SB4 star. The radial velocity variations of its four components can be explained by two separate 
eclipsing binaries. In contradiction to previous photometric findings, we show that the observed composite spectrum as well as the 
derived masses of all four of its components correspond to spectral type F. 

Conclusions. The observed small dip in the light curve is not caused by a transit-like phenomenon but by the eclipses of the second 
binary system. We find evidence that KIC 4247791 might belong to the very rare hierarchical SB4 systems with two eclipsing binaries. 
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1. Introduction 

The Kepler mission was designed and launched with the primary 
aim of searching for transiting pla nets, and indeed, ha s found 
hundreds of planetary candidates (iBorucki et alJ |20I Th and a 
large number of eclipsing binaries (EBs). In addition to those 
already very interesting objects, eve n more unusual syst ems 
such as multiple trans iting planets (jL issauer et al. 20 1jJ) or 
triply eclipsing trip les (ICarter et al.l [201k .Derekas et aLj l201ll: 
ISlawson et al.ll20I ih have been discovered. 

In their analys is of low-mass EBs i n the Kepler QO and 
Ql data releases, ICoughlin"et hereafter Cll, find 

that one object does not belong to the main-sequence, KIC 
4247791 = 2MASS J 19083956-h3922369. The authors determ- 
ine its temperature to be reff=4 063K, an orbital period of 
4^100866, a combined mass of 1.28 M©, and a combined radius 
of 3.82 Rq. According to these values, the system must contain 
at least one evolved component. Moreover, a periodic transit-like 
feature is superimposed on the EB light curve. With respect to 
the main eclipse, this feature drifts with time (Fig. [T]) and the 
authors derive possible periods of either 2^02484 or 4^04969. 
According to CI 1, the feature could be explained in the case of 
the shorter period by (1) a background EB with no visible sec- 
ondary, (2) a circumbinary transiting object, or (3) a transiting 
object around one of the stars in an almost 2:1 resonant orbit 
with the binary. In the case of the longer period, the explanation 
could be (4) a fore- or background EB with 0.98752 times the 



* Based on observations with the 2-m Alfred- Jensch-Telescope of the 
Thiiringer Landessternwarte Tautenburg. 



main period of KIC 4247791 that has nearly identical depths of 
the primary and secondary eclipses. 

KIC 4247 791 also has an e ntry in the Kepler eclipsing bin- 
ary catalogue (iPrsa et al.ll201 1|) and both EBs are included in its 
second release (Isiawson et al.ll201 ll) . 

We obtained 15 high-resolution spectra covering the orbital 
period of KIC 4247791 well and use them in this investigation 
to find the real cause of the observed transit-like feature in the 
Kepler light curve of the star and to derive tighter constraints 
on the fundamental parameters of its components. For this pur- 
pose, we measure the radial velocities (RVs) of the components 
of KIC 4247791, analyse the averaged composite spectrum, de- 
rive the orbits of the found binary systems, and combine the res- 
ults with those obtained from the Kepler light curve analysis. 



2. Observations 

All spectra have been taken with the Coude-Echelle spec- 
trograph attached to the 2-m Alfred- Jensch-Telescope of the 
Thiiringer Landessternwarte Tautenburg. The spectra cover the 
wavelength range 4 700 A to 7 400 A with a resolution of 32 000. 
The exposure time was 40 min per spectrum. The dates of ob- 
servation are listed in Table lA.ll Spectrum reduction was done 
using standard ESO-MIDAS packages including bias and stray- 
light subtraction, filtering of cosmic rays, flat-fielding using a 
halogen lamp, optimal order extraction, wavelength calibration 
using a ThAr lamp, and normalization to the local continuum. 
Nightly shifts in the instrumental zero point were corrected by 
using a large number of telluric O2 lines. 
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Figure 1. Detrended Kepler light curve (Q0-Q3 data) phase fol- 
ded to 4^? 100871. A shallow dip drifts with time. Around 55 
cycles are covered, each shifted upwards by 2 mmag (times the 
cycle number) for visibility. 
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Figure 2. Observed spectrum in the region (black, BJD= 
2455 819) and best-fit synthetic spectrum (red) computed for 
Teff =6 400 K, log g=3.7, [M/H]=-0.2 dex, and v sin /=40 km s'^ . 



3. Spectrum analysis 

In a first step, we determined the predominant spectral 
type from one observed, composite spectrum showing the 
sharpest (unsplit) lines. For the determination of the stel- 
lar p^ameters^we used a method identical to that described 
in iLehmann et al.l (1201 ll) . to derive values of reff=6400K, 
log g = 3.1, V sin / = 40 km s"^ and a metallicity of -0.2 dex re- 
lative to the solar value. Figure[2l compares an observed spec- 
trum with the best-fitting synthetic spectrum for the 11J3 region. 
The results are representative of the composite spectrum. This 
spectrum has to be decomposed to derive the spectral types of 
the single components. The derived temperature corresponds to 
spectral type F and contradicts the value of 4 063 K given by 
CI 1. In particular, the observed B aimer line profiles are incom- 
patible with such a low temperature value. 

4. Orbital solutions from cross-correlation 

Cross-correlation functions (CCFs) were obtained from the 
cross-correlation of the observed spectra with the synthetic tem- 
plate spectrum derived in Sect.[3]but using v sin / = 1 km s"^ in- 
stead of 40 km s"^ to obtain the maximum spectral resolution in 
the CCFs. In most of the CCFs, we discovered four compon- 
ents. The positions of these components represent the RVs of 
four stars in two binary systems with orbital periods close to 
the periods of 4^100866 (the two stronger peaks in the CCFs) 
and 4^049687 (the two smaller peaks) given in CI 1. This can be 
seen clearly if we arrange the CCFs according to the positions in 
orbital phase based on the two periods (Fig. [5]). 
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Figure 3. CCFs vertically shifted with orbital phase based on 
periods of 4^100866 (top) and 4^49687 (bottom). The peak po- 
sitions belonging to the corresponding orbit are marked by aster- 
isks, and the positions of minimum and maximum RV separation 
by plus signs. 
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Figure 4. The Fe 1 4 957 A line observed at two different orbital 
phases showing strong (continuous line) and almost no (dashed 
line) splitting. 

The line splitting can also be seen in the single spectra. 
FigurelH compares the Fel 4 957 A line from two spectra taken 
at different orbital phases where the noise was reduced by av- 
eraging over 9 wavelength bins. One spectrum clearly shows a 
splitting into four components belonging to system I, its stars A 
and B, and system II, its stars A and B, whereas the other spec- 
trum has almost no splitting. 

The RVs were determined from the CCFs using multiple 
Gaussian fits. From 14 spectra, we could determine the RVs of 
all components. We used four Gaussians in 12 cases and five 
Gaussians in the 2 cases where we detected a fifth contribution 
produced by weak telluric lines. In one spectrum, we could only 
resolve the RVs of the brighter EB. Table lA.ll lists the RV values 
obtained for the four components. 

The orbital solutions were calculated usi ng the metho d of 
differential corrections to the orbital elements (ISterndl 19411) . We 
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Figure 5. From top to bottom: RVs of components I A (circles) 
and IB (crosses) folded with the period of 4^100871, 0-C re- 
siduals of I A and IB, RVs of components II A (circles) and IIB 
(crosses) folded with the period of 4^049732, and 0-C residuals 
of II A and IIB. 

assumed circular orbits in each case. For both systems, we fixed 
the orbital periods and times of y- velocity passage/primary mi- 
nimum (iterated with those obtained in Sect.O. Figure[5] shows 
the RVs of both systems and the corresponding 0-C values fol- 
ded with the two diff'erent periods. Three outliers can be seen 
for component lA, one for IB and one for IIB. These outliers 
were rejected when calculating the orbital solutions. They can 
arise when the RV components coincide and are poorly resolved 
in the CCFs, in particular at/around phases and 0.5, which 
moreover can be aff'ected by the Rossiter eff'ect (for A and B, 
respectively). The derived orbital elements are listed in the up- 
per part of Table[T] 

Phase zero in Fig. [5] corresponds to the time of primary mi- 
nimum (rMiN in Table[T]). We see that in both systems the RVs 
of the more massive components A that have slightly smaller 
RV amplitudes, decrease at this phase. This means that the more 
massive components are eclipsed at Turn and should have higher 
surface brightnesses. It is remarkable that not only the periods of 
the two systems are almost identical (but significantly diff'erent 
as can also be seen from their Kepler light curves) but that also 
the RV semi-amplitudes K and the mass ratios q are very similar. 

5. Light curve analysis 

The Kepler database provides publicly available data for 228 
days (2009-05-02 to 2009-12-16; quarters Q0-Q3) of photo- 
metric observations of KIC 4247791 (Fig. 1), including 88 days 
taken in short cadence mo de (SC, Q2, 1 niin expo sure time). 
On the basis of JKTEBOP feouth worth et al.ll2004allbh . the fight 
curve was analysed simultaneously for both EBs using the IDL 
implementation mpf it of the Levenberg-Marquardt algorithm. 
A fit was made to the two binary light curves using 14 free para- 



Table 1. Parameters and spectral types (SpT) derived from the 
spectroscopic, light curve, and combined analysis. Turn gives 
BJD-2 455 000, errors in units of the last digits are given in par- 
entheses. 





system I 


system II 




A B 


A B 


K[km s-'] 


91.9(4) 101.6(3) 


88.6(3) 91.4(4) 


y [kms~^] 


-16.2(9) -17.1(8) 


-11.8(2.0) -11.7(1.7) 




0.905(7) 


0.970(7) 


P[d] 


4.100871(4) 


4.049732(71) 


Tmin [BID] 


1.1453(1) 


1.353(1) 


/ [deg] 


79.5(1) 


80.0(8) 


ecoso; 


0.00012(6) 


0.0020(3) 




-0.00570(4) 


0.0009(4) 


e 


0.00570(4) 


0.0022(3) 


CO [deg] 


178.7(6) 


65(24) 


R/a 


0.1570(4) 0.1503(4) 


0.100(9) 0.092(8) 




0.576(1) 0.424(1) 


0.55(1) 0.45(1) 


Lab I ^i+ii 


0.76(2) 


0.24(1) 


a[AU] 


0.0742(2) 


0.0681(2) 




2.50(1) 2.40(1) 


1.46(13) 1.35(12) 


Mab [Mo] 


3.24(3) 


2.56(3) 


M[Mo] 


1.70(2) 1.54(1) 


1.30(2) 1.26(2) 


logg [dex] 


3.87(1) 3.87(1) 


4.22(8) 4.28(8) 


SpT 


FOIV F2IV 


F7V F8V 



meters in each case: binary period, phase zero-point, ratio and 
sum of the radii, inclination, surface brightness ratio, eccentri- 
city (^sino;, ^cos6l)), quadratic limb-darkening, and reflection 
efl'ect coefficients. The mass ratio was fixed to the values determ- 
ined from the RV measurements. The relative brightness ratio as 
well as a contribution from a (constant) fifth background light 
was also taken into account. 

The Kepler light curves are influenced by instrumental ef- 
fects, which were treated separately for each quarter or even for 
sections separated by an interruption of the observations. This 
is fitted simultaneously with the binaries by piecewise second 
or third order polynomials. Starting from an initial guess for a 
total of 52 parameters using identical stellar parameters, zero 
eccentricity, zero fifth light, and neglecting systematic trends, 
we employed the IDL routine mpf it to perform an optimisa- 
tion, mpf it calls JKTEBOP twice to evaluate the contribution 
of each binary. Since several of the parameters are closely re- 
lated, the result of the optimisation depends partly on the initial 
conditions. After a 5-(T clipping to remove outliers, we there- 
fore started 50 optimisation runs where the initial conditions 
were randomly varied around the first result. The one with the 
lowest ;^2 was taken as the final result, mpf it also provides 
errors from the diagonal of the co variance matrix. For more 
realistic errors in the presen ce of correlated no ise, we how- 
ever, followed the approach of lSouthworthI (l2008l) . To take non- 
Gaussian residuals into account, a series of 500 synthetic light 
curves was constructed using the best-fit model adding arbit- 
rarily shifted residuals as noise. Since the residuals from the 
long and short cadence data differed significantly, they were 
treated separately. These light curves were again fitted using 
various start values. From the nearly Gaussian distribution of 
the parameters around the best-fit relation, we determined the 
1-(T uncertainty limits. These are typically an order of magnitude 
larger than the formal errors provided by mpf it 

The derived parameters are listed in Tabled and Fig. [6] illus- 
trates the quality of the fit. In addition, we verified the results 
by fitting each EB separately. For this purpose, the data were 
decomposed iteratively by phase folding the light curve to the 
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Figure 6. Disentangled and phase-folded light curves. Left: 
Light curve (EB II model subtracted) fitted by the EB I model. 
Right: Light curve (EB I model subtracted) fitted by the EB II 
model. The corresponding 0-C residuals are shown in the lower 
panels. 

periods of the EBs and calculating the medians in one-minute 
bins. Many cycles had to be be covered to minimise the residual 
correlation eff'ects introduced by the ellipsoidal variations and 
the similar periods of both EBs. The Q2-SC data covers about 
21 orbit cycles of EBI and 27% of the full beat cycle (324.5 
days) of both EBs. 

With the known periods, inclinations, and RV amplitudes of 
each component, we were able to derive the stellar radii and 
masses listed in Table[T] From the masses and radii, we rou ghly 
estimated the spectral types using the tables of [Gmvl(l2005l) . We 
derived a spectral type F for all four stars, which agrees with 
the observed composite spectrum. For main-sequence stars, we 
would expect radii of about 1.5 R© for the stars of system I and 
of 1.3 R© for the st ars of system II and logg of about 4.3 for 
all the components (lGrayll2065b . Comparing these values with 
Table[Tl we see that both components of system I must be slightly 
evolved. 

From the derived orbital period and the stellar radii, we were 
able to estimate the v sin / of the components of system I assum- 
ing synchronized rotation. We measured 31 kms"^ for compon- 
ent A and 30 km s"^ for B. This is compatible with the 40 km s"^ 
derived from the composite spectrum in Sect.[3] that we should 
recall is aff'ected by an additional broadening since this spectrum 
was taken close to but not exactly during the (double) y- velocity 
passage. 

6. Discussion 

Our results clearly show that the eff'ects observed in the Kepler 
light curve can be explained by two EBs and not by some transit- 
like phenomenon. It is difl&cult to understand, however, that all 
the similarities in the spectroscopically derived parameters of 
the two binary systems are caused by chance and one has to 
ask wether there is a possible connection between the systems. 
Therefore, the most interesting question is whether both EBs are 
gravitationally bound and form a quadruple system. 

To exclude visible blends, we did a careful search for any 
contaminants in the Kepler pixel data and in available ground- 
based images of the corresponding field of view. KIC 4247791 
can be found in one J-band image of the UKIRT InfraRed 
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Figure?. Left: 1.12'x 1.50'section of a Kepler full-frame image 
(FFI). The used extraction aperture for KIC 4247791 is marked 
by the red borders. Right: UKIRT /-band image of the same re- 
gion. The Kepler aperture is overlaid. The marked, nearby blend 
is KIC 4247807. 

Deep Sky Survey (UKIDSSQ, DR WSERV4, see lLawrence et al.l 
r2007), as well as in the Sloan Digital Sky Survey (SDSS-Ilfl 
DR8) in all five (u, g, r, /, z) bands. The UKIRT image has the 
higher spatial resolution with a pixel scale of O'.'2/pixel but the 
image of KIC 4247791 itself is overexposed. The SDSS provides 
data of a'396/pixel and KIC 4247791 is non-saturated in u and 
z, where the Gunn-z image has the better seeing. We used the 
UKIRT /-band image to check for possible blends in the vicin- 
ity of the star and the SDSS z-band image to investigate its flux 
profile. 

An inspection of the /-band image and a comparison with 
the pixel data of the Kepler images identified an object (KIC 
4247807, Kepler magnitude = 14.808 mag) separated by 
1472 from KIC 4247791 (Kp = 11. 260 mag) close to the aper- 
ture used for extraction (Fig .IT]). However, the following finding 
excludes the possibility that KIC 4247807 is EBII: the spectra 
were taken with a 2'' slit, where the slit orientation on the sky 
varies with the hour angle of the telescope. For any object that 
did not lie completely inside the slit, we were unable to derive 
as precise a RV curve as we obtained from the spectral features 
of EB II. The close object is 3.5 mag fainter than KIC 4247791. 
Assuming that 10% of the fight of KIC 4247807 from the PSF 
wings falls into the aperture, its fifth light contribution should be 
at a level well below 1%. 

Figure[8] shows in its upper left panel a contour plot of the z- 
band image of KIC 4247791 . To enhance the resolution, we con- 
structed the mean point spread function (PSF) from all brighter 
and non- saturated star images in the z-band field by applying a 
shift and add, sampling all images to ten times higher resolution 
than provided by the original ones. Two-dimensional Gaussian 
fits where used to eliminate all subimages showing any asym- 
metry in the flux profiles. We ended up with a PSF built from 63 
star images with a 1708 FWHM, the same value as we obtained 
from the image of KIC 4247791 itself. The upper right panel in 
Fig. [8] shows the image of the star de convolved w ith the PSF 
using the Lucy-Richardson algorithm (iLucyll 19741) . It is almost 
perfectly symmetric, possessing a FWHM of 0^.^32. For compar- 
ison, we simulated a nearby star of 0.2 times the brightness of 
KIC 4247791 at a distance of Off 5 using the derived PSF. The 

^ http://surveys.roe.ac.uk/wsa 
^ http://www.sdss3.org 
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Figures. Top: Contour plots of the SDSS z-band image of 
KIC 4247791 (left) and of the image deconvolved with the PSF 
(right). The position scale is in arcsec. Bottom: The same with a 
simulated faint companion in (7/5 distance. 



result is shown in the lower panels of Fig.[8j The deconvolved 
image is clearly asymmetric having a FWHM of Off 44 in the dir- 
ection towards the faint companion and of Off 35 in the perpen- 
dicular direction. We conclude that the separation between EB I 
and EB II must be closer than Off 5, which still implies that they 
may be either affected by an accidental blend or be gravitation- 
ally bound systems. 

There are some findings that are consistent with a quadruple 
configuration of KIC 4247791 such as the luminosity ratio that 
provides similar distances of both EBs or the similar inclinations 
and masses that are indicative of a common formation. The sim- 
ilarity between the systemic RVs also implies that it is a phys- 
ical associated system. The small diff'erence in the system RVs 
of --5 km s"^ could, however, be due to the orbital motion of both 
EBs around each other and might be used to derive an upper limit 
to the separation between both EBs. For a circular and edge-on 
orbit, this upper limit would be 1200 years. We note that this 
limit is proportional to ^K~^ meaning that if we observe e.g. 
only a quarter of the maximum RV diff'erence the period would 
shorten to about 20 years. In the latter case, the corresponding 
variations in the systemic y- velocities of the two systems should 
be safely detectable from observations spanning about two years. 



7. Conclusions 

Using the period derived in CI 1 from the minima in the Kepler 
light curve, we have been able to model the RV variations of 
the two main components identified in the CCFs derived from 
our newly obtained spectra. Two additional features have been 
observed in almost all the CCFs, which vary with the period of 
4^049751 that was proposed by Cll as a possible period to ex- 
plain the occurrence of a shallow dip in the Kepler light curve 
drifting with time against the main minima. 

Our spectroscopic analysis has uncovered four stars mov- 
ing in two separate binary systems. The basic evidence for such 
an interpretation comes from the time-dependent splitting of the 



spectral lines into four components that we have been able to as- 
sign in a unique way to four stellar components moving in two 
separate orbits owing to the small diff'erence between the two 
periods being significant, and visible as small moving dips in 
the Kepler light curve. 

Our results for the combined photometric and spectroscopic 
observations are consistent with the scenario of two binaries, 
both lying along the same line of sight, showing eclipses, having 
almost the same orbital periods, RV amplitudes, and mass ratios. 
This seems to be a very rare if not unlikely scenario and raises 
the questions of whether these two systems are physically con- 
nected. The observed similarities may be indicative of a quad- 
ruple system, although we cannot entirely rule out a moderate 
back- or foreground scenario. 

All derived fundamental parameters of EB I are well defined. 
For the primary and secondary, respectively, we obtained masses 
of 1.70+0.02 Mo and 1.54+0.01 M©, and radii of 2.50+0.01 R© 
and 2.40+0.01 Ro and logg of 3.87+0.01 for both compon- 
ents, corresponding to slightly evolved stars of early F-type. 
For EBII, the derived masses have similar accuracy to that 
of EB I but the errors in the radii and log g are distinctly lar- 
ger. We obtained Ma =1. 30+0.02 M©, = 1.46+0.1 3 R©, and 
log gA =4.22+0.08 for the primary and Mg^l. 26+0.02 M©, 
7?fi=1.35+0.12Ro, and loggfi=4.28+0.08Mo for the secondary 
corresponding to late F-type main sequence stars. 

To our knowledge, only two stars have been found so far 
to host two EBs in hierarchical systems, the visual binaries 
ADS 9537 AB and V994 Her. ADS 9537 AB was discovered by 
iBatten & Hardid (1 1965b to consist of two WUMa-type EBs 
where the two sys tems orbit e ach other with a period of about 
22000yr (iBatten & Lu 19861). iLee et"aD (l2008h discovered that 
the SB4 star V994 Her hosts two EBs consisting of B8+A0 and 
A2+A4 main sequence stars. From astrometric data, they estim- 
ate a period of the visual binary of few thousand years. T here 
are a few other remarkable objects. IShkolnik et al.l (l2008l) dis- 
covered a SB4 quadruple system, BD-22 5866, consisting of one 
(K7+K7) EB and a second (M1+M2), non-EB system, orbiting 
each other with a period shorter than 10 yr. OGLE 05134 3.14- 



691837.1 was proposed as a SB4 candidate \ yith two EBs bvlOfii 



(2008) b ut has since been f ound to be SB 2 (iKolaczkowski et al 



.2010: Riviniusetal.l 1201 Ih. The SB 3 quadruple system V379 
Cep was found by iHarmanec et all (l2007l) to host two binaries 
where one is an EB. 

Further observations would be necessary to confirm the 
quadruple nature of KIC 4247791, namely data extending the 
present time base and designed to detect the orbital motion 
between systems I and II based on transit timing or RV vari- 
ations. In a forthcoming paper, we intend to decompos e the 
spectr a of the components using the KOREL program ( Hadraval 
l2Q06h and we will analyse in greater detail the decomposed light 
curves and spectra to provide more stringent constraints on the 
fundamental stellar parameters, in particular the eff'ective tem- 
peratures of the components, and hence the evolutionary state 
and distance modulus of KIC 4247791. 
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Appendix A: Measured radial velocities 



Table A.l. RVs in km s ^ derived from the four components in 
the CCFs where aRV gives the mean measurement error. 



BJD-2455 000 


RViA 


RVm 


RViiA 


RVim 


700.529717 


7.82 


-46.25 


60.25 


-85.15 


702.571156 


-39.50 


10.92 


-86.85 


62.50 


725.484044 


40.41 


-93.13 


66.07 


-95.22 


726.433284 


53.99 


-93.06 


-36.43 


10.44 


734.477064 


65.83 


-108.09 


-28.07 


5.81 


754.407123 


67.44 


-110.81 


13.18 


-37.77 


757.429167 


-60.46 


30.93 


73.98 


-98.91 


782.488513 


2.61 


-37.39 


49.85 


-74.03 


787.374417 


81.94 


-117.03 


-60.90 


52.13 


788.514326 


-15.77 


-15.77 


-75.66 


53.74 


793.387013 


-102.52 


76.62 


29.65 


-53.81 


794.419142 


-49.49 


19.83 


66.29 


-90.44 


796.454312 


19.35 


-55.60 


-87.92 


68.29 


796.483258 


15.31 


-51.26 


-86.30 


66.50 


819.341840 


-12.25 


-12.25 






ARV 


0.15 


0.15 


0.28 


0.24 



6 



